Anharmonic effects on infrared spectra of GaAs and GaP 
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Abstract. With a study of bulk GaAs and GaP, we report detailed infrared spectra and their temperature dependence from 
first principles. Fine features and low-temperature trends in the predicted response are confirmed by experiment. The spectra 
are calculated from the far-infrared through twice the frequency of the zone-center optical phonon. The treatment relies on 
recent developments in the theories of dielectric polarization and the phonon-phonon interaction. 



Ab initio prediction of the optical response of simple 
solids over a very broad spectral range is a significant 
achievement of theory and is founded on a quantitative 
understanding of the electronic structure Ql. In a related 
advance, complete characterizations of lattice dynamics 
in the harmonic approximation, where phonons do not 
interact, are now routine |2]. 

The harmonic approximation, however, is insufficient 
for a close inspection of the frequency dependence of the 
dielectric response in the infrared (IR). For instance, in 
the absence of direct multiphonon transitions, the har- 
monic approximation predicts infinitely sharp absorp- 
tion at OrO' the frequency of the long-wavelength, trans- 
verse optical phonon, or the dispersion oscillator, and 
nowhere else |3l0|. Measured spectra, however, show fi- 
nite linewidth and asymmetric lineshape of the IR-active 
phonon and very broad background absorption, which 
can exhibit fine features and singularities ™d 
temperature-dependent absorption in the far-infrared 01 • 
These effects can be explained by incorporating the 
phonon-phonon interaction into the theory of the optical 
response, where anharmonic terms in a crystal's interi- 
onic potential cause the phonons to interact 
As an example of this, an infrared photon can, through 
anharmonic corrections to the IR-active phonon's quan- 
tum state, be absorbed by a pair of phonons whose mo- 
menta sum to zero and whose energies sum to that of 
the photon. These sum processes endow the spectra with 
rich structure. At nonzero temperatures, the onset of dif- 
ference processes, where one phonon is emitted and an- 
other is absorbed along with the photon, give a strong 
temperature dependence to the far-infrared absorption. 

Realistic, observable quantities arising from phonon- 
phonon interactions are now calculated as well as mea- 
sured, such as the linewidth, asymmetric shape and tem- 
perature dependence of Raman modes L12iil3i.il4i.il5i.il6t 



fl?!. In this Letter, the treatment of the phonon-phonon 
interactions, or anharmonic effects, closely follows the 
theoretical studies of Raman spectra, where knowledge 
of a phonon self-energy is of central importance. Exam- 
ination of the phonon self-energy's effects in the dielec- 
tric response is interesting in its own right, and the abil- 
ity to calculate it may aid in the design of infrared op- 
tical components. Phonon-induced dipoles are given by 
the Born effective charges. These charges enhance the 
longitudinal optical phonon frequency while coupling 
the transverse optical phonon to infrared photons, and 
they are calculated with the modern theory of polariza- 
tion its', 1?', . With these charges and the IR-active 
phonon's self-energy, the spectral profile and the oscil- 
lator strength can be calculated, and all the optical con- 
stants can be obtained. 

Figures 1 and 2 summarize our results and their 
comparison with measured spectra. Their examination 
suggests that two-phonon processes are present over 
the entire infrared spectral range through 2(Oto- Below 
coto, we predict the measured |8] low-energy absorp- 
tion tail well into the far-infrared spectrum, and the ex- 
pected temperature dependence. The lineshape of the IR- 
active phonon is strikingly reproduced. Between cOjo 
and 207-0, absorption and transparency features demon- 
strate a convincing correspondence with measured data. 
These are associated, respectively, with critical points 
and gaps in the zero-momentum, two-phonon density 
of states. There is more structure to the two-phonon 
spectrum of GaP than of GaAs, including narrow bands 
of suppressed absorption, because its larger ionic mass 
mismatch widens the energy difference between various 
phonon branches. The data for GaAs are from Ref. |6] 
and references therein. Some data points are from fits 
to an oscillator model, and are "generally not deter- 
mined better than ± 30 %" [6]. The reststrahlen and low- 
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FIGURE 1. The one-phonon (line) and zero-momentum 
two-phonon (+) densities of states are plotted in (a) for 
GaAs. In (b) and (c) the calculated absorption coefficient 
is plotted at K and 300 K, respectively. Plotted theoreti- 
cal spectral features have been moved to 0.98 of their cal- 
culated position, in cm^', to facilitate comparison with the 
experimental 300 K spectra, which are plotted in (c) (+). 
The arrows index, from left to right, the two-phonon critical 
points: TA(L) + TA{L), TA{X) + TA(X), TO{L) + TO{L), and 
LO(r) +LO(r). The optical response is characterized with the 
absorption coefficient, cc{v) = 4zvk{v), where k{v) is the in- 
dex of absorption, and v is photon energy in cm^'. The data 
are from Ref. [6]. 



FIGURE 2. Same as in Fig. 1, for GaP, with the exception 
of the overtone at L, as noted in the text. The data are from 
Ref. 1 5] and Ref. 1 7] (x) as described in the text. 



and their many-body treatment, and finally present re- 
sults and details of our calculations. 

In GaAs and GaP, infrared optical properties are dom- 
inated by phonons. The Hamiltonian density for a polar 
crystal in the presence of an electric field, E, can be writ- 
ten as 



frequency data for GaP are from Ref. 10] and references 
therein. The low-temperature and room-temperature data 
points above the reststrahlen are from Ref. 15], and the 
data are from transmission measurements on a high- 
purity sample. 

The extent of similarity between predicted and ob- 
served spectra strongly indicates that modern first- 
principles calculations and perturbative formalism reflect 
the microscopic processes governing infrared optical re- 
sponse, and may serve to stimulate further work toward 
detailed characterization of crystalline materials. 

Below we summarize the role of phonons in the di- 
electric response, discuss the anharmonic contributions 



(1) 

where and co^ are the normal coordinate and fre- 
quency, respectively, of a phonon with wave vector and 
branch indexed by X. The unit cell volume is Q. and the 
number of Bravais lattice sites is A^. The dimensions of 
the quantity P'""' are dipole moment per unit volume. 
Anharmonic potential energy terms are assumed to be 
small, and expressed with //„„/,. Direct dipole transitions 
to multiphonon states are represented with H,„„. This 
coupling has been observed and calculated L31.. ,32.1 . but 
is not accounted for here. 



The macroscopic polarization in Eq. Q can be written 



as. 



v=TO Tj 



(2) 



In the above, is the mass of ion T, the phonon nor- 
mal coordinate index, v, runs strictly over the IR-active 
phonons, and e^j is a component of the polarization of 
such a phonon. Cartesian directions are indicated with 
the indices / and /'. The charges are expressed as the 
change in macroscopic polarization with respect to sub- 
lattice displacement; 



3 pmac 



(3) 



Modern theory gives elegant and tractable prescriptions 



with which to evaluate Zjj [18, 19]. 

Formally, the following expressions can be derived 



with a formalism analogous to that of Kubo 12 ir . A 
heuristic oscillator model for the response can be writ- 
ten by expressing the effects of //„„/, as a complex damp- 
ing term, or self-energy that modifies (Oro- ^ro(f^) = 
A(a)) — iyico). Substituting Eq. (|3 into Eq. {I), solving 
the driven-oscillator equation of motion for an IR-active 
phonon's normal coordinate, and then taking the sum in- 
dicated in Eq. the polarization is found to be: 



(4) 



v=TO 



1 



CO + COto + '^To{(o) co ~ coto -'^to{o}) 



The dielectric function follows immediately. It is writ- 
ten ESQ: 



v=TO 



1 



1 



CO + COto + '^Toio)) co- coto-'^to{o3) 



(5) 



with e^y representing the high-frequency dielectric con- 
stant. When all the IR-active modes are symmetry- 
related, as is the case in a zincblende system, their fre- 
quencies and self-energies are equal, and denoted by the 
subscript TO above. 

Low-temperature phonon-phonon interactions have 
been the subject of thorough field-theoretic analyses 0, 
[13 El HimiH. Modern implementations suc- 
cessfully calculate realistic thermodynamic and spectral 
properties which originate from such interactions lEsl 
129*1. demonstrating the valuable predictive nature of first- 
principles studies. Phonon-phonon scattering events are 
ranked with a dimensionless expansion parameter, for 




FIGURE 3. Lowest-order processes contributing to the IR- 
active phonon self-energy. The vertex of diagram (a) represents 
a four-phonon interaction and is related to energy shifts associ- 
ated with quartic terms in lattice strain. Diagram (b) expresses 
shifts associated with thermal expansion. Only diagram (c) con- 
tributes frequency-dependence to the self-energy. 



which Van Hove has posited the ratio of zero-point ion 
delocalization to nearest-neighbor lengths |24]. The the- 
ory establishes lowest-order contributions to the anhar- 
monic self-energy of a long-wave phonon from three 
distinct multiphonon processes l24l |25L E6ll . These are 
shown in Figure 3. One of these processes, (a), involves 
one four-phonon vertex (second order in the expansion 
parameter), and two, (b) and (c), involve three-phonon 
vertices (each first order in the expansion parameter) and 
are second-order in //„„/,. The four-phonon process and 
one of the three-phonon processes represent simultane- 
ous absorption and reemission of a phonon by the IR- 
active phonon. Such interactions are instantaneous, and 
therefore independent of frequency, contributing exclu- 
sively a line shift 1 30]. Our interest is in spectral structure 
and multiphonon processes, so we have neglected these 
frequency-independent terms. In any case, our calcula- 
tions assume the theoretically relaxed rather than actual 
lattice constant, so the positions of cOto related spec- 
tral features are affected by other factors as well. 

The remaining term, (c), reflects the IR-active 
phonon's coupling to phonon-pair states of opposite 
crystal momenta, k and — k, where their role has been 
compared to that of a "viscous medium" |3], and the 
matrix element between the IR-active phonon and such 
a state is: 



Mafik^ E 



dDV, 



du 



(6) 



The derivative of the dynamical matrix, (k), is taken 
with respect to the sublattice relative coordinate. Phonon 
branches are denoted a and j3 . For a zincblende crystal, 
the relevant anharmonic contribution to the Hamiltonian 



is: 



E E ^cePk?ka?-k/3?v. (7) 



The prefactor includes the reduced mass, = (w^go 
ni^l)^"^, or n — {mqI^ +mp^)^^. The anharmonic inter- 
action in Eq. ([7) leads to a frequency-dependent shift and 



damping of the IR-active phonon, expressed as the setf- 
energy: 



-TO 



l+«ka + n-k/3 \ / 1 



(8) 



(Oka 



"k/J - 



(Oka 



(O-kp 
nka - «k/3 



co + irj 



^(Oka - (Okp - (0 - It] J XCOkp - COka - (0 - 11] 

which is derived in detail in Refs. . The sum- 

mation includes mode frequencies, (Oka, and Bose oc- 
cupation factors, Hka- We work within a frozen-phonon 
framework |34, 35], using density-functional theory, 
pseudopotentials and a plane-wave basis I3a.l37ll38ll3^ . 
The results for the two distorted crystals are subtracted to 
evaluate the derivative in the dynamical matrix of Eq. |7} 
at 122 points within the Brillioun zone. Fitting 22 in- 
equivalent real-space interatomic force constants facil- 
itated increasing the Brillouin zone sampling to ^ 10*' 
points. 

In Figures 1 and 2, along with the theoretical and ex- 
perimental values for the absorption coefficient, we have 
also plotted the density of momentum-conserving two- 
phonon states for comparison with the absorption fea- 
tures. The dispersion-oscillator-two-phonon quasiparti- 
cle model we have adopted is shown to agree remarkably 
with experimental results near the reststrahlen and qual- 
itatively away from it. The arrows in the figures repre- 
sent a few sum-process critical points as identified in the 
caption. The wave vectors are indicated with the stan- 
dard notation for the high-symmetry points of a face- 
centered cubic reciprocal lattice, and the branches are la- 
beled L(T)A(0) for longitudinal (transverse) and acous- 
tic (optical). The three most energetic modes are denoted 
optical, the doubly degenerate modes denoted transverse, 
and the singlets denoted longitudinal. The identified two- 
phonon overtones are identical for each material, with the 
exception of the optical overtone at L, where the singlet 
is more energetic in GaP, and the reverse is true for GaAs. 
The doublet overtone is indicated for GaAs, and the sin- 
glet for GaP. The longitudinal-transverse band crossing 
in GaAs is a result of its similar ionic masses, as dis- 
cussed below. Features above 20^0 are evident in the 
theoretical spectra of GaP because of slight overbending 
in the calculated phonon dispersion. 

Above the reststrahlen, the theory reproduces spectral 
structure better than it does overall amplitude. For in- 
stance, the calculated absorption for GaAs is typically 
between one half and one quarter of the measured value 
between 400 cm^' and 600 cm^^ The relative values 
for the measured absorption are more reliable than their 
absolute values |6|, and this fact may contribute to the 
discrepancy. But it may be that the two-phonon fea- 



tures are not completely described through the IR-active 
phonon self energy, and that direct multiphonon transi- 
tions should be considered. As can be seen from Eq. (jSJl, 
our model cannot give absorption above Icojo, where 
three-phonon and higher-order processes must be con- 
sidered [33]. 

The temperature dependence of the absorption below 
100 cm^' is ascribable to difference processes, where 
a phonon is absorbed and one of the same momentum 
but differing branch is emitted. Such processes are ex- 
pressed in the third and fourth terms on the right side 
of Eq. (|8}, and are only possible at temperatures corre- 
sponding to finite phonon occupation numbers. The ab- 
sorption grows as the material heats and the occupation 
numbers increase, and this behavior has been measured 
in GaAs for three far-infrared frequencies ] 8 ] . 

The GaP spectrum possesses a richer two-phonon 
structure than that of GaAs, including narrow bands of 
weak absorption, a consequence of its larger ionic mass 
mismatch, which widens the energy difference between 
various phonon branches. In the case of similar masses, 
such as those in GaAs, the optical and acoustic branches 
approach near degeneracy at high-symmetry wave vec- 
tors, thus narrowing the gap between the bands associ- 
ated with the branches. The same is true for the trans- 
verse and longitudinal optical branches, which, in the 
case of GaAs, show greater dispersion and again nearly 
degenerate behavior at high-symmetry points. These ef- 
fects can be seen through comparison of phonon densi- 
ties of states of Figure 1 and Figure 2, as well as through 
comparison of the absorption. Figure 2 shows that the 
primary transparency band we calculate for GaP, just 
below 700 cm^\ is observed. Although this trough in 
the absorption appears weak next to the sharp IR-active 
phonon feature, it constitutes a local variation in the ab- 
sorption by two orders of magnitude. For both materi- 
als, spectral structure is also evident in the temperature- 
dependent far-infrared absorption below the reststrahlen, 
where critical points can be seen in the theoretical plots. 

In short, we have obtained quite realistic results for the 
infrared spectra and their temperature dependence from 
first principles with an anharmonic two-phonon theory. 
The behavior above twice the IR-active phonon cannot 
be modeled without higher-order multiphonon processes. 
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